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Abstract 

Within the cell, biopolymers are often situated in constrained, fluid environments, e.g., cy- 
toskeletal networks, stretched DNAs in chromatin. It is of paramount importance to understand 
c/5 1 quantitatively how they, utilizing their flexibility, optimally respond to a minute signal, which is, 

in general, temporally fluctuating far away from equilibrium. To this end, we analytically study 
4h . viscoelastic response and associated stochastic resonance (SR) in a stretched single semi-flexible 



chain to an oscillatory force or electric field. Including hydrodynamic interactions between chain 



segments, we evaluate dynamics of the polymer extension in coherent response to the force or field. 

in '. 

We find power amplification factor of the response at a noise-strength (temperature) can attain 

Ov 

the maximum that grows as the chain length increases, indicative of an entropic stochastic reso- 

OO " 

nance (ESR). In particular for a charged chain under an electric field, we find that the maximum 
also occurs at an optimal chain length, a new feature of ESR. The hydrodynamic interaction is 
found to enhance the power amplification, representing unique polymer cooperativity which the 

X 

fluid background imparts despite its overdamping nature. For the slow oscillatory force, the res- 
onance behavior is explained by the chain undulation of the longest wavelength. This novel ESR 
phenomenon suggests how a biopolymer self-organizes in an overdamping environment, utilizing 
its flexibility and thermal fluctuations. 
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I. INTRODUCTION 



Within the cell, the biopolymers such as DNA, proteins, and actin filaments are confined 
subject to various kinds of forces or fields. In chromosome, for example, double stranded 
(ds) DNA is condensed in a fascinating hierarchy of varying structures, interwound via his- 
tone proteins subject to tension in chromatin level Associated with a cellular membrane, 
cytoskeletal network consists of actin filaments, upon which myosins bind to exert tension. 
An universal feature that underlies in these systems over various levels is the polymer flexi- 
bility that can adapt to ambient thermal noises as well as the external forces. Another one 
is the aqueous background which transmits hydrodynamic interactions (HI). In this vein, we 
pose a fundamental question : How the nature allows the confined (stretched) biopolymer 
systems to utilize the ambient thermal fluctuations and the system flexibility for maximum 
cooperativity to an external signal in such an overdamping media? 

The stochastic resonance (SR) is a noise-induced phenomenon, where the background 
noise can enhance coherence and resonance of a non-linear system to a small additive periodic 
signal 

a a. 

The SR is a measure of coherence, synchrony, and sensitivity of the system 
to external influences. The exemplar model is a bistable system of a single particle that 
undergoes interwell hopping transition coherent to the signal at an optimal noise-strength 

as 

. Also, the SR is recently found to occur even in a non-linear monostable system of 



single particle or single degree of freedom 



As linearly interconnected systems they are, polymers manifest unusual cooperativity and 
collective dynamics in response to those external stimuli. It has recently been shown that 
the chain interconnectivity amplifies the thermally-induced hopping and the associated SR 
of a flexible polymer moving under a bistable potential js, 7|. The SR for such polymeric 
systems is of entropic nature, which shows the maximal coherence to an oscillatory force 
not only at optimal noise-strengths, but also, remarkably, at the optimal chain lengths that 
depend on chain flexibility and conformational changes. This entropic SR suggests novel 
possibilities for manipulating the dynamics of single biopolymer transitions utilizing their 
flexibility. 

To address the SR of biopolymers in constrained situations mentioned earlier, in this work 
we study the dynamics of an end-tethered worm-like chain (WLC) subject to a constant 
tension and an additional small oscillating force. The end-to-end distance (extension) of 
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FIG. 1: The coordinates of the semi-flexible chain in a fluid. 

the chain is indeed found to respond most coherently at an optimal noise-strength, and, at 
optimal length for the case of the polyelectrolyte under an electric field. It turns out that 
the most pronounced SR is due to the longest-wave length lateral undulation of the chain 
in response to the oscillation. We show and discuss how these system parameters and HI in 
water with temperature-dependent viscosity affect the SR behavior. Eventually, the effective 
dynamics for the chain extension is governed by the Langevin equation subject to a non- 
linear monostable free energy due to the constant tension force and entropic restoring force, 
as well as a non-constant drag coefficient that incorporates HI and chain inextensibility. 

Sec. II is a brief introduction of the dynamics of a stretched semi- flexible chain in linear 
response to a time dependent force, including HI. In Sec. Ill, we specialize in the response of 
the chain extension, the SR in particular, to a small oscillating force applied to a chain end. 
SR is evaluated as a function of the noise-strength as well as chain parameters and HI. In 
Sec. IV, we extend the analytical results of Sec. Ill to the case of a charged polymer under 
an oscillating electric field. In Sec.V, to analyze the nature of entropic SR and effect of HI, 
the chain dynamics is projected into one-dimensional dynamics of the extension under a free 
energy and an effective drag. 



II. DYNAMICS OF SEMI-FLEXIBLE CHAIN UNDER A TIME-DEPENDENT 
TENSION 

Consider a semi-flexible chain whose one end is pulled by a time-dependent force. Pro- 
vided that the force is strong enough to stretch the chain nearly to the full contour length L, 
we represent the position vector of a segment at an arc length s in a cylindrical coordinates 
r(s,t) = (h(s, t), z(s, t)) where h(s,t) is radial undulation vector and z(s,t) is arc length 
projected along the force (Fig. [p. 

Due to segment inextensibility, |^f^| = 1, and small undulation, (^^) < 1, we 
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consider here, the chain extension along the force is 
Z(t)~- 



ds ( |^ 

OS 



l- 1 - 

2 



L , fdh(s) 

as 



ds 



(1) 

'0 \ us J * Jo \ / 

Of our major interest is how the average end-to-end distance, (Z(t)), responds to a time- 
dependent force denoted by /(£), depending on the chain parameters (such as stiffness and 
length) and on the fluid environment characteristics (such as the temperature and HI). In 
this study we confine ourselves to the force f(t) = fo + Sf(t) where \Sf(t) \ <C f . Initially(at 
t = 0), the chain is brought to a thermal equilibrium under a constant and uniform force 
/o- The dynamic response of the mean chain extension to a time-dependent force in a 
viscous liquid was addressed in references si, 9]. For logical flow we recapitulate here the 
basic equations and obvious expressions that were already derived or known. However, a 
modification is made in the chain mobility by incorporating HI as well as self-frictions. This 
is crucial for assessing the effect of HI on SR in our study. 

To study the dynamics of chain extension afterwards, we consider an effective Hamiltonian 
for the chain pulled by the force given by 



2 /„ dSK [ ds" 



ds 



-f(t)z(L,t) 
d 2 h(s" N 2 



ds 2 



+ f(t) 



dh(s] 
ds 



-f(t) (L-z(0,t)) 



(2) 



where two terms in the [ ] represent chain bending and stretching energy, and the last term 
is trivial with z(0,t) =0. k is the bending rigidity of the chain and is related to the chain 
persistence length L p by k = ksTLp. The dynamics of the h(s,£) is given by the Langevin 
equation 

dh(s,t) rL 
in 



Here, A(s,s') 



U(s 



Of 
s') + 



[ ds'A(s,s') ■ F(s',t) + R(s,t). 
Jo 



(3) 



i 



87"?|r as /| 



1 + 



\*ss'\ 2 



is the mobility tensor consisting of self 



friction (with the friction coefficient ( per unit length of the chain) and hydrodynamic 
interaction between two different segments s and s', at a distance r ss > = r(s) — r(s'), where 
1 is unit tensor, 77 is the fluid viscosity. The R(s, t) is the Gaussian and white noise. F(s' , t) 
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is the density of the radial force acting on segment s', 



F(s',t) = -6H/5h(s',t) 



ds'* J J w \ds' 2 

Within the small undulation approximation we consider, the A(s, s') = 
1 [6(s - s')/( + l/8n7]\s - s'\) to reduce the Eq. as 

By the Fourier transform, e.g., h(q, t) = ds e iqs h(s,t), the Eq. §5§ is rewritten {§]] as 

^M = - W (g,t)h(g,t) + R(g,t), (6) 

where oo(q,t) = A(q) [nq A + f(t)q 2 ], A(q) = l/£ + K (qa)/4irr), and K (qa) is the modified 
Bessel function of the second kind with a cutoff a introduced as the chain diameter 10]. 
The average chain extension is given by 

q=Tr/a 



(z(t)) 



l ~VT2 E <? 2 <h(^)-h(-^)>- (7) 



L 2L 2 

q=ir/L 



From Eq. (Q, we obtain to the order linear in Sf, 

t 

(h(q, t) ■ h(-q, t)) = (h(q, 0) • h(-q, 0)> - J dt'm(q, t - t')5f(t') 

o 

where (h(q, 0) ■ h(-q, 0)) = 2k B TL/(Kq A + f q 2 ) and 



t-t 



Ak B TLA(q) 
m(q, t-t)= e . 

[Kq 2 + fo) 

The 

r q = [2u(q, 0)]- 1 = 1/ [2A(g) {nq 4 + f q 2 )] (8) 

is the relaxation time of undulation mode q stretched by the force fo. At t — 0, Eq. (JTj) 
yields the well-known relation between equilibrium stretching force and extension Zq/L = 
(Z(0))/L = 1 - k B T/^f or 

'-SO"*)" W 
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This is due to entropy change associated with the extension. The Eq. (J7J) is rewritten as 

(6Z(t)) ee (Z(t)) - Z 
t 

= J dt'M(t-t')Sf{t'), (10) 



o 



in terms of a dynamic response function, 

M{t-t') = ^2q 2 m(q,t-t') 



1 r' a 2 4k B TLA(q) 
2vr J n / L [nq 2 + f ) 



(11) 



The Eq. ( flOj) and Eq. ( flTT) represent the chain extension in linear response to time- 
dependent tension of general form in a viscous fluid under HI. We mention here that these 
equations were used in relaxation dynamics of the chain extension under time-dependent 
drag exerted by an end-anchored large bead [8| and also in viscoelastic dynamics under 



oscillatory force without HI [9|, respectively. Our study here differs from theirs in the explicit 
expressions for the mobility A(q) = 1/C + Ko(qa) / 'Airq , and also in the approximation via a 
linear response theory described next. 

Although we will analyze Eq. ( flTT) in detail for SR, we need to obtain its approximation 
that allows us to gain much simpler analytical understanding. We first note that the Eq. fflUj) 
is the consequence of a dynamic linear response theory, according to which, 

M(t) = -^ T j t (6Z(t)5Z(0)} , (12) 

where (...)o means the average over the equilibrium ensemble in the absence of Sf(t). 
Integrating the above relation, we obtain the time correlation function 

(8Z(t)5Z(0)) = -k B T [ dt'M(q,t~t') + ((5Z) 2 } 

Jo 



(k B T) 2 L r ,a , 1 

dq- — -^e T i . 



J^L (nq 2 + f o y 



(13) 



which at t — is reduced to equilibrium fluctuation 

{k B T) 2 L r /a 

TT Lil 1 {nq 2 + f 



((SZ)\ = ^-L^ I dq - o _ x2 . (14) 




4 6 8 10 
t [10" 6 sec] 



FIG. 2: The correlation function of the chain extension, (5Z(t)5Z (0))q (solid curve). We use 
T = 37°C, L = 100 nm, L p = 50 nm, a = 2 nm and /o = 0.5 pN. The single relaxation time 
approximation (dZ(t)8Z(0))o ~ ((5Z) 2 )oe~ t / T is depicted by a dashed curve. 



Fig. [2] depicts the relaxation behavior of the correlation function (Eq. (ITS]) ) of a chain 
stretched by the force fo for the parameters typical in stretched DNA. The correlation decays 
faster in time as the stretching force or the bending rigidity increase. In the long time regime, 
the correlation function is well approximated as 

(SZ(t)SZ(0)) ~ ((5Z) 2 ) e-^, (15) 

where r is the longest relaxation time associated with the mode of a longest wavelength 
(q = ir/L) undulation given by 

T= 2A(7r/L)((7r/L)4« + (7r/L)2/o)' ^ 
This single relaxation time approximation is valid in the long time limit as well as at t = 0, so 
that it can incorporate exactly the response to the static and low-frequency forces, although 
it yields a rather slowly decaying correlation (Fig. [2]). In this approximation M(t) is given 
as 

M{t) = ^^-e-^, (17) 
k B l t 

in a form much simpler than Eq. (TTTT) . 
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III. THE ENTROPIC STOCHASTIC RESONANCE OF CHAIN EXTENSION 
UNDER AN OSCILLATING FORCE ON AN END 



Consider a time-dependent force in the form, 5f(t) = f m sin(Qt), to be applied at one end 
of the chain while the other end is held fixed. Since we assume that the tension is uniformly 
acting along the chain, we should consider the cases where the tension propagation time, 



Tf = r]L 2 y/Kl '(25L p f^ 2 ) [l2| should be much shorter than the time scales of the f(t) and 
the (5Z(t)). From the Eq. f lTUj) . the long time asymptotic behavior of the chain extension 
becomes 

(5Z(t)) = f m [x' sin(ta) - X " cos(fif)] , (18) 

where 



„ Ak B TL r' a nr q q 4 A 2 (q) 

x - _ l aq 2 + Q2 ■ 



71 



tt/L lq 



(19) 



The Eq. ( !T8l) is rewritten as 

(6Z(t)) = \x(mfmSm(nt-<p), (20) 



where \x\ = a/ x' 2 + X" 2 is the absolute magnitude of dynamic susceptibility x(^) an d 
tp = tan -1 (x" /x') is the phase delay. The power amplification factor, P = \x\ 2 , is a measure 
of coherence or synchrony with phase delay to the oscillatory force. 

Although we will study the above P for SR behavior, its single relaxation time ap- 
proximation P is useful. Noting that the dynamic susceptibility x(^) is nothing but the 
Fourier-Laplace transform of M(t), x(fi) — Jo°° dt e int M(t), 

((5Z) 2 ) 1 

X{n) = k B T (21) 
which, in the limit ft — > 0, recovers the static susceptibility of the chain xo — ((dZ) 2 ) /k B T 
in response to the static force fo- We then obtain the power amplification factor in this 
approximation, 

r f((«y>°y 1 (22) 

F \ k B r ) i + (Tfir {U] 
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(b) L = 118nm, L = 50 nm 
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FIG. 3: The power amplification factor P vs. the temperature in Celsius for three frequencies 
= 10 s _1 (solid), $1 = 10 4 s" 1 (dashed) and Q, = 10 5 s _1 (dotted), with the hydrodynamic 
interactions (HI) included. For the semi-flexible chain with dsDNA parameters, a = 2 nm, L p = 50 
nm, and Jq = 0.5 pN, the P are maximized at optimal temperatures indicated by the arrows, for 
contour lengths (a)L = 100 nm, (b)L = 118 nm, (c)L = 122 nm, and (d)L = 125 nm. 



We examine characteristics of the power amplification factor P for a semi-flexible chain 
taking the parameter values of a double stranded DNA (dsDNA) fragments, a = 2 nm, 
L p = 50 nm and 20 < L < 125 nm. In addition, we incorporate temperature dependence 
of the water viscosity 77 given by the empirical relation log = ai ( 20 ~ T ^~°2 ( 20 ~ T ) where 
r]o = 1 cP is the viscosity at T = 20°C, a x = 1.1709, a 2 = 0.001827, a 3 = 89.93 and T is 
temperature in the Celsius scale [3]. A chain with L = 100 nm and L p = 50 nm pulled by 
fo = 0.5 pN in water at T = 310 K = 37°C is stretched by 0.86L < (Z) < 0.96L, according 
to the Eq. Q). Fig. [3] depicts the P(T) with four different contour lengths L = 100, 118, 
122 and 125 nm, respectively. The angular frequencies Q of the driving force are chosen to 
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be smaller than the tension propagation frequency Qf = 2n/Tf, which is about ~ 10 7 s _1 
for the parameters we consider. 

As indicated by arrows in the Fig. EJ the P{T) for a semi-flexible chain (with parameters 
considered) the dsDNA can have peaks with respect to the noise-strength T, which is the 
main characteristics of the SR, for the certain ranges of polymer parameter L and angular 
frequency Q. In one- dimensional picture of the chain dynamics which will be described in 
Sec.V, the chain extension is subject to an effective non-linear potential or the free energy, 
which depends explicitly on temperature due to the entropy of the chain. Consequently, 
the SR obtained here is of the entropic nature and is dubbed as the Entropic Stochastic 
Resonance (ESR). 

The P(T) behavior in the Fig. Ovaries sensitively to small variations of L. For example, 
the solid curves (the P(T) with the lowest angular frequency Q = 10 s _1 ) in the Fig.^a), (b), 
(c) and (d) are drastically distinct from one another only with difference of L by ~ 10 nm; 
they monotonically decreasing, non-monotonically varying with a peak, and monotonically 
increasing, respectively, within the range of temperature shown in the Fig. |3j In addition, 
each P(T) in the Fig. |3^a), (b), (c) and (d) changes greatly as Q decreases from 10 5 s _1 
to 10 4 s _1 , while it negligibly changes as Q decreases below 10 4 s _1 , converging rapidly to 
the Q — > limit, which is explained below and seen in Fig. Ufa) later. We find that the 
response of a 118 nm-long chain with dsDNA parameters pulled by 0.5 pN force and driven 
by a minute oscillation at Q = 10 4 s _1 is maximized at an optimal noise-strength which is 
around the physiological temperature, 37°C (dashed curve in the Fig. [31(b)) . Of particular 
note is that the SR peak is quite broad (rather than feeble) within the temperature range 
where the water can exist in the liquid state, suggesting that a biopolymer in such condition 
can be widely responsive to a minute external force. 

If Q <C r _1 , the power amplification factors, coached by Eq. (1221) . become independent of 
dynamics details and are largely determined by the static susceptibility xo = ((^^) 2 )o/^b7 1 , 
from which we can evaluate the conditions for SR peaks. Performing the integration Eq. (IT^l) 
for the cases where foL 2 / (ksTLp) 3> 1, 




(23) 
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FIG. 4: Enhancement of the power amplification P by HI. (a) P vs. O with (solid curve) and 
without HI (dashed curve), at T = 37° C and L = 100 nm. HI can enhance the amplification about 
twice at Q = 10 6 s _1 (see dotted vertical lines). The power amplification P by the single relaxation 
time approximation is depicted by the dash-dotted curve, (b) P vs. T with and without HI, at 
Q = 10 6 s" 1 and L = 100 nm. (c) P vs. L/L p with and without HI, at T = 37°C and O = 10 6 s -1 . 
Each cross point of the dotted lines depicts the P at fi = 10 6 s _1 , T = 37°C and L = 100 nm. All 
the other parameters are the same as those used in the Fig. [3l 



leading to the optimal temperature for SR, 

^ 9 f L 2 



(24) 



256 ksLp 

Among others, this shows that T* increases (sharply) with L, in consistency with the peaks 
at low frequencies in Fig. [3](b) and (c), as mentioned earlier. It establishes an important 
condition for SR; the chain does not have a SR if fo = 0, i.e., unless the chain is constrained. 
Also, it indicates that neither an ideal chain (L p = 0) nor a rigid chain (L p = oo) can 
support the SR. 

Fig. H] shows enhancement of the power amplification P(Q, T, L/L p ) by HI. The Fig. Ufa) 
depicts P vs. Q with and without HI, at T = 37°C and L = 100 nm. With HI turned off, 
A(g) — > 1/Q with ( put as 37T7], the P is reduced to about half that at Q = 10 6 s -1 , meaning 
that HI enhances the P. Note that the effect of HI is not significant for the low angular 
frequencies below r _1 , which is about 10 5 s _1 for the parameters considered over which the 
chain segmental dynamics associated with HI is rapidly relaxed out and the static response 
(f2 — > 0) dominates. This feature is explained by the single relaxation time approximation 
P, which is depicted by the dash-dotted curve in FigH^a). It evidently become an excellent 
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approximation to P as f2 decreases to zero. Since the relaxation time r decreases with HI 
as shown by Eq. ffTB"]) . the Eq. f l22|) indicates SR is enhanced by HI for Q > r _1 . Since HI 
makes r _1 larger by the factor ~ K^ijia/ L) f=s | ln(7ra/L)| compared with that without HI, 
HI effect can be greatly enhanced for longer contour length, as shown by Fig. Wis). The P 
at fixed Q = 10 6 s -1 , T = 37°C and L = 100 nm with and without HI are depicted as cross 
points of the dotted lines, in the Fig. HJ 



IV. THE ENTROPIC STOCHASTIC RESONANCE OF CHAIN EXTENSION 
UNDER AN OSCILLATING ELECTRIC FIELD 

In this section, we consider a polyelectrolyte chain which is uniformly charged with a 
linear charge density A, subject to a time dependent electric field E(t) = — [E + <5E(t)} z 



with 5E <C Eq. The effective Hamiltonian is then 

k ( d 2 h(s 



U 



ds 



ds* 



-E(t)\Z(s,t) 



(25) 



An effect of electrostatic interactions between the segments are incorporated into a renormal- 
ization of the persistence length [l^, l3] that depends upon the salinity of the background 
fluid. The perturbation effect of the 5E(t) on the Hamiltonian is 



W = -LX5E(t)Z cm (t), 



(26) 



where Z cm (t) = J Q L dsZ(s,t)/ L is the chain center of mass position along the field. Hence, 
unlike the previous case of an end-stretching force, the Z cm (t) is the coordinate conjugates to 
the force -L\5E(t) acting on entire segments. The mean of Z cm (t) in the small undulation 
limit is shown to be half the mean end-to-end distance. According to phenomenological 
and analytical arguments, equilibrium end-to-end distance under the strong electric field is 
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15 



16] 




(27) 

which means that the uniform electric force corresponds to the end-stretching force via the 
relation f Q = L\X\E Q /A. The response of the mean center of mass position to the oscillating 
electric field 5E{t) = 5E sin (Qt) is then obtained by replacing f m by L\X\5E /A and x by x/2 
in the Eq. ( |T9l and Eq. ( 1201) . The power amplification factor then becomes Pg = |x| 2 /4 with 
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FIG. 5: The power amplification factor Pe vs. temperature in Celsius for three frequencies f2 = 10 
s _1 (solid), VL = 10 2 s _1 (dashed) and O = 10 3 s _1 (dotted), with HI included. For the stiff chain 
taking parameters of a F-actin filament with a = 6 nm, L p = 15 fim and A = —1 e/nm, the Pe can 
have maxima at optimal temperatures indicated by the arrows for contour lengths, (a)L = 1.5 //m, 
(b)L = 2.2 /im, (c)L = 2.3 /mi and (d)L = 4.1 /mi. 



fo replaced by L\\\E / 4. Note that the electric field induces a tension which is proportional 
to the chain length L. 

Under an electric field E Q = 50 V/cm, a tethered stiff chain with parameters of F-actin 
filament fragment, L > 1 fim, L p = 15 /im, A = — 1 e/nm, and a = 6 nm, is extended 
as much as Z /L > 0.9 following the Eq. (1271) . Considering a small additive driving field 
5E sin (Qt) , we plot the -Pe(T) for the cases of the stiff chain extension at four different 
contour lengths, in Fig. [51 including the hydrodynamic interactions. This also shows that 
the ESR occurs for optimal temperature, similar to the case driven by the stretching force. 
In the slow driving regime (solid curves in the Fig. [5](a), (b), (c) and (d)), the Pe also show 
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FIG. 6: The Pe vs. L/L p with the various O at T = 37°C, with HI included. The Pe has the peak 
with respect to a contour length, which becomes sharper for smaller f2. All the other parameters 
are the same as those used in the Fig. [5j 

distinct behaviors within the temperature considered, which are monotonically decreasing, 
non-monotonically varying with a peak, and monotonically increasing, respectively. 

Using the single relaxation time approximation P E = (((SZ) 2 ) /2k B T) 2 /(l + (rfi) 2 ). At 
low frequencies, temperature peaks are due to by the static susceptibilities. Following the 
argument behind Eq. ([24"]) . the optimal temperature is obtained via the replacement fo 
by L\X\E /A in Eq. (j24l). as T* = (9/256) \X\E L 3 /{4k B L p ), which is consistent with the 
dependence of T* on L, in Fig. |5]for low frequencies. 

Remarkably, as depicted in Fig. |6] for a fixed T, there exists a sharp peak of the power 
amplification Pe(L/L p ) at an optimal contour length, which is a new feature of the ESR. This 
is due to the competition between the undulational flexibility and the electric-field induced 
stretching rigidity. As the chain flexibility increases with contour length, the Pe increases 
as in the end-tethered case, while it decreases due to enhanced chain rigidity induced by the 
electric field effect on the chain segments. As the angular driving frequency Q decreases, the 
ESR (magnitude of the peak) sharply increases. Note that, in the Fig. El the optimal contour 
length (for the peak) is around 4.1 fim, seemingly independent of the angular frequencies. 
The replacement of fo by L|A|£^o/4 in Eq. (123]) yields the optimal contour length for the SR 
peak 




(28) 
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FIG. 7: Enhancement of the power amplification Pe by HI. (a) Pe vs. f2 with (solid curve) and 
without HI (dashed curve), at T = 37°C and L = 2.2 /mi. HI can enhance the amplification by 
almost twice at O = I0 3 s _1 (see dotted vertical lines) (b) Pe vs. T with and without HI, at 
Q = 10 3 s- 1 and L = 2.2 /mi. (c) P E vs. L/L p with and without HI, at T = 37°C and Q = 10 3 
s _1 . Each cross point of the dotted lines depicts the Pe at Q = 10 3 s _1 , T = 37°C and L = 2.2 /mi. 
All the other parameters are the same as those used in the Fig. [5j 



For the parameters considered therein, L* w 4 /im in agreement with the Fig. |6] for the 
r _1 . For this optimal length, the Fig. |5^d) shows Pe{T) increases monotonically with 

T. 

In Fig. [7J we also show HI effect on the Pe for the field-driven case, where its dependency 
on Q (or T) is similar to the tension-driven case. We find that the magnitude of Pe{Q) at 
T = 37°C and L = 2.2 fim can be enhanced by almost twice that at Q > r _1 ~ 10 2 s _1 , 
due to HI, as depicted in Fig. [7(a). Fig. [TJb) and Fig. [7(c) depict Pe vs. T and P vs. 
with and without HI, respectively for the parameters indicated. Remarkably, as shown in 
the Fig. [TJc), the ESR peak can be highly magnified by HI at the optimal length. 



V. THE EFFECTIVE 1-D DYNAMICS OF SEMI-FLEXIBLE CHAIN EXTEN- 
SION AND ENTROPIC STOCHASTIC RESONANCE 

Conventionally, theoretical study of SR has mostly considered one (or few) dimensional 
stochastic system, for example, a Brownian particle under a one-dimensional (ID) double 
well potential. On the other hand, hitherto, we have investigated the SR of a stretched 
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chain by considering the motion of many segments, which are locally governed by the elastic 
interactions under condition of the segmental inextensibility as well as the hydrodynamic 
interaction. To understand entropic and dissipative nature of SR, we construct the effective 
1-D dynamics of the chain extension Z(t), based on the single relaxation time approximation. 

First the Eq. (JH]) allows us to derive an effective potential or free energy of the chain 
extended by the constant tension f , 

k B T L 



HZ) 



4L p 1 - Z/L 
k B T ( L 
IE 



-foZ 



l- z/l (i-z Q /Lyj- 



(29) 



Here Z = (Z(0)) is the equilibrium extension where the free energy is minimum. The ^(Z) 
as depicted in Fig. [HI is anharmonic, rising rapidly as the chain is fully extended, due to the 
chain entropy (the first term in Eq. (1291)). 

The 1-D chain dynamics we consider is essentially that of the Brownian motion fluctuating 
around the free energy minimum. Using the approximation of a single relaxation time, we 
now establish the effective dynamical equation of the slow chain extension. It is assumed to 
be the Langevin equation of the form, 

d7{Z) 



TZ 



dZ 



+ 5f(t)+£(Z,t). 



(30) 



Here T is the hydrodynamic drag coefficient of the whole chain under extension Z. It is 
related to the random force £,{Z,t) via a fluctuation-dissipation theorem, 



{£(Z,t)Z(Z,Q)) = 2k B Tr6(t). 



(31) 



Since, in our approximation, (Z(t)) relaxes exponentially toward Z Q with the characteristic 
time r in the absence of Sf(t), we have 



r 0J T{Z) 



(32) 



Matching the above with Eq. ( 1301) allows us to find the Z-dependent drag, 



V{Z) 



Z-L(l 



k B T ■ 



knT 



t k B T 



Z-Z 4L F 



[AL p (l-Z/Lf 
1 1 
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;i - Z/Lf (1 - Z /Lf 



(33) 
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FIG. 8: The free energy J-(Z) (solid curve) and effective friction coefficient T{Z) divided by 
7 = 67T7/(a/2) (dashed curve) as functions of the extension Z. The dotted vertical line depicts Zq 
(~ 80 nm) where free energy is the minimum. We use T = 37° C, L = 100 nm, L p = 50 nm, a = 2 
nm and /o = 0.5 pN. 



In the Fig. El T(Z)/ r y is depicted by the dashed curve with segmental friction 7 put as 
6in](a/2). T(Z) rapidly increases with Z above Z w 80 nm as depicted by the dotted 
line. This is mostly due to effect of the chain entropy that prohibits the chain to extend to 
the full length, in this effective model. It should be noted also that HI manifests through 
r, Eq. (I16p . by making it, consequently the chain drag (Eq. ( 133]) ) smaller by the factor 
~ K^ajV) rs I ln(7ra/L)| compared with that without HI. Consequently, the chain exten- 
sion dynamics, effectively reduced to a single degree of freedom Z, is stochastic subject to a 
anharmonic potential (non-linear force) and a non-constant friction due to the chain entropy 
and hydrodynamic interactions of the undulating, extended chain. The extension oscillates 
around Z in maximum coherence to an oscillating force at SR condition. Compared with 
the SR conventionally studied in ID bistable systems with a constant friction, there are two 
points worthy of remarks: First, the monostable free energy arising the chain connectiv- 
ity and entropy manifest the SR (ESR). Second, another feature of the cooperativity, the 
hydrodynamic interaction, can enhance it by reducing the drag. 
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VI. SUMMARY AND DISCUSSION 



We analytically have studied dynamical response and associated stochastic resonance 
(SR) in a stretched semi-flexible chain to an oscillatory (AC) field or force. Incorporating 
the chain segmental elasticity and inextensibility, we set up an effective Hamiltonian of 
the chain in terms of the chain transverse undulation field which is assumed to be small, 
and evaluated the undulation and extension in response to the external oscillation. From 
the dynamic susceptibility of the chain extension, we directly obtained power amplification 
factor as a measure of the SR. The chain extension is then found to be most synchronous and 
coherent to the oscillation at an optimal noise-strength (temperature), revealing an SR of 
entropic nature. In particular, for the case of an applied AC field we found a new type of the 
entropic SR (ESR), where the power amplification is maximum at an optimal chain length 
due to competition between associated conformational entropy (chain flexibility) and electric 
field (chain rigidity). This ESR is enhanced by the hydrodynamic interactions between the 
chain segments, due to its cooperativity even in overdamping media. 

On the other hand, from the way the chain extension responds to a time- dependent field 
or force, we have constructed the time correlation function of the stretched chain extension 
in its absence, using the linear response theory. It decays more rapidly than the exponential 
decay with the slowest relaxation of time associated with the transverse undulation of longest 
wave-length. Yet the latter describes quite accurately the long time behavior of the stretched 
chain dynamics and the SR to low frequency oscillation. We have constructed the relevant 
1-D dynamics of the stretched chain extension, which is found to be subject to a monostable 
free energy due to entropic force and a drag reduced by the hydrodynamic interaction, both 
of which depend upon the chain extension. This means that stretched polymer is a non-linear 
monostable system to support a novel kind of ESR. The simple 1-D stochastic equation we 
derived (Eq. (132))) should be useful to simulate the stretched semi-flexible chain dynamics 
in correlation with experiments. 

The biopolymers in vivo are usually subject not only to thermal noise (temperature) but 
also to athermal, non-equilibrium noises of various kind due to active cellular environments. 
A non-equilibrium noise is usually characterized by a broadband frequencies, but as we 
have shown, the component corresponding to the lowest frequency is most amplified at 
an optimal temperature or an optimal contour length. Although our semi-flexible chain 
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model may be too simple to represent the real biopolymers in action, our result on SR may 
suggest biological systems and their environments have evolved to adopt the physiological 
temperatures for optimal biological functions, as recently have been shown in the SR of 
ion-channel The novel features of ESR for a single biopolymer we discussed here can 
)e tested by single molecule experiments for stretched biopolymers such as tethered dsDNA 



20[ and stiff actin-flaments 21] under oscillating fields. 
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